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ABSTRACT 

Infrared  Space  Observatory  Long  Wavelength  Spectrometer  grating  observations  of  the  ultraluminous  infrared 
galaxy  Arp  220  shows  absorption  in  molecular  lines  of  OH,  H20,  CH,  NH,  and  NH3,  as  well  as  in  the  [O  i]  63  //m 
line  and  emission  in  the  [C  n]  158  pm  line.  We  have  modeled  the  continuum  and  the  emission/absorption  of  all 
observed  features  by  means  of  a  nonlocal  radiative  transfer  code.  The  continuum  from  25  to  1300  pm  is  modeled 
as  a  warm  (106  K)  nuclear  region  that  is  optically  thick  in  the  far-infrared,  attenuated  by  an  extended  region  (2") 
that  is  heated  mainly  through  absorption  of  nuclear  infrared  radiation.  The  molecular  absorption  in  the  nuclear 
region  is  characterized  by  high  excitation  due  to  the  high-infrared  radiation  density.  The  OH  column  densities  are 
high  toward  the  nucleus  (2-6  x  1017  cm-2)  and  the  extended  region  (~2  x  1017  cm-2).  The  H20  column  density  is 
also  high  toward  the  nucleus  (2-10x  1017  cm-2)  and  lower  in  the  extended  region.  The  column  densities  in  a  halo 
that  accounts  for  the  absorption  in  the  lowest  lying  lines  are  similar  to  what  are  found  in  the  diffuse  clouds 
toward  the  star-forming  regions  in  the  Sgr  B2  molecular  cloud  complex  near  the  Galactic  center.  Most  notable 
are  the  high  column  densities  found  for  NH  and  NH3  toward  the  nucleus,  with  values  of  ~1.5  x  1016  cm  2  and 
~3  x  1016  cm-2,  respectively,  whereas  the  NH2  column  density  is  lower  than  ~2  x  1015  cm-2.  A  combination  of 
photodissociation  regions  (PDRs)  in  the  extended  region  and  hot  cores  with  enhanced  H20  photodissociation 
and  a  possible  shock  contribution  in  the  nuclei  may  explain  the  relative  column  densities  of  OH  and  H20, 
whereas  the  nitrogen  chemistry  may  be  strongly  affected  by  cosmic-ray  ionization.  The  [C  n]  158  pm  line  is  well 
reproduced  by  our  models  and  its  “deficit”  relative  to  the  C  n/far-IR  ratio  in  normal  and  starburst  galaxies  is 
suggested  to  be  mainly  a  consequence  of  the  dominant  non-PDR  component  of  far-infrared  radiation,  although 
our  models  alone  cannot  rule  out  extinction  effects  in  the  nuclei. 

Subject  headings:  galaxies:  abundances  —  galaxies:  individual  (Arp  220)  —  galaxies:  ISM  — 
galaxies:  starburst  —  infrared:  galaxies  —  radiative  transfer 


1.  INTRODUCTION 

With  a  redshift  of  z  =  0.018,  Arp  220  (IC  4553/4)  is  the 
nearest  and  one  of  the  best-studied  ultraluminous  infrared 
galaxies  (ULIRGs).  The  tails  observed  in  the  optical,  together 
with  the  double  highly  obscured  and  compact  (0"3)  nuclei 
observed  in  the  near-  and  mid-infrared,  as  well  as  in  the  mil¬ 
limeter,  strongly  suggest  that  the  enormous  luminosity  of 
Arp  220,  ~1012  Lq,  is  the  result  of  galactic  merging.  Never¬ 
theless,  the  concrete  physical  process  responsible  is  still  a 
matter  of  debate:  the  proposed  sources  are  hidden  active  nuclei 
and/or  bursts  of  star  formation. 


1  Based  on  observations  with  the  Infrared  Space  Observatory,  an  ESA 
project  with  instruments  funded  by  ESA  member  states  (especially  the 
principal  investigator  countries:  France,  Germany,  Netherlands,  and  the 
United  Kingdom)  and  with  the  participation  of  ISAS  and  NASA. 

2  Visiting  Astronomer,  Harvard-Smithsonian  Center  for  Astrophysics,  60 
Garden  Street,  Cambridge,  MA  02138. 

3  CSIC,  IEM,  Departamento  Astrofisica  Molecular  e  Infrarroja,  Serrano 

123,  E-28006  Madrid,  Spain. 


Molecular  observations  of  Arp  220  provide  unique  clues  to 
the  physical  and  chemical  processes  occurring  in  the  nuclei 
and  their  surroundings.  In  the  millimeter  region,  CO  obser¬ 
vations  have  been  carried  out  with  increasingly  high  angular 
resolution  (Radford  et  al.  1991a;  Scoville  et  al.  1991,  1997; 
Okumura  et  al.  1994;  Sakamoto  et  al.  1999).  In  particular  they 
have  shown  that,  on  the  one  hand,  CO  emission  arises  from  a 
region  significantly  more  extended  than  the  nuclei  (~3"-4"), 
and  on  the  other  hand  that  the  CO  (2-1)  to  CO  (1-0)  intensity 
ratio  is  lower  than  1,  thus  suggesting  that  CO  mainly  traces 
low-density  regions  (<103  cm-3).  Observations  of  molecules 
with  high  dipole  moment,  like  CS  and  HCN,  have  revealed 
that  the  fraction  of  molecular  gas  contained  in  dense  clouds 
[»(  h2]  >  104  cm  3)  is  much  larger  than  in  normal  galaxies, 
yielding  ~1010  M0  of  dense  gas  (Solomon  et  al.1990,  1992). 
Radford  et  al.  (1991b)  found  that  HCN  (1-0)  and  HCO+  (1-0) 
peak  strongly  toward  the  nuclei  but  also  show  low-level 
extended  emission.  More  recently,  Aalto  et  al.  (2002)  have 
detected  emission  from  the  high-density  tracers  HNC  and 
CN,  and  the  relatively  low  HCN/CN  and  HCN/HNC  intensity 
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ratios  were  attributed  to  widespread  photodissociation  region 
(PDR)  chemistry. 

The  launch  of  the  Infrared  Space  Observatory  (ISO)  opened 
a  new  window  for  the  study  of  the  physical  and  chemical 
properties  of  ultraluminous  infrared  galaxies.  Despite  the 
lack  of  angular  and  spectral  resolution,  the  observations  of 
Arp  220 ’s  far-infrared  spectrum  from  40  to  200  pm  (Fischer 
et  al.  1997,  1999)  and  of  a  number  of  individual  lines  in 
the  Short- Wavelength  Spectrometer  range  (Sturm  et  al.  1996) 
provided  new  insights  in  our  understanding  of  the  ionic,  atomic, 
and  molecular  content  of  the  galaxy.  These  wavelength  re¬ 
gions  are  of  great  interest,  because  the  bulk  of  the  enormous 
luminosity  is  emitted  in  the  far-infrared,  and  also  because 
they  contain  lines  of  important  molecular,  ionic,  and  atomic 
species.  Skinner  et  al.  (1997)  reported  the  detection  of  the 
35  pm  OFl  line  in  Arp  220.  Fischer  et  al.  (1997,  1999)  found 
that  the  far-infrared  molecular  absorption  lines  of  OH, 
H20,  CH,  and  NH3  are  significantly  stronger  in  Arp  220  than 
in  less  luminous  infrared-bright  galaxies,  while  the  fine  struc¬ 
ture  lines  from  ionic  species  are,  to  the  contrary,  extremely 
weak.  Luhman  et  al.  (1998,  2003)  found  that  in  Arp,  like  in 
other  ULIRGs,  the  [C  n]  to  far-infrared  luminosity  ratio  is 
typically  nearly  an  order  of  magnitude  weaker  than  in  lower 
luminosity  infrared-bright  galaxies.  Sturm  et  al.  (1996)  re¬ 
ported  the  detection  of  two  ortho-H2  pure  rotational  lines,  in¬ 
dicating  that  high  masses  of  gas  are  subject  to  PDR  conditions 
and/or  shock  activity.  Although  the  few  other  ULIRGs  for 
which  either  full  Long  Wavelength  Spectrometer  (LWS) 
spectra  or  partial  scans  in  the  vicinity  of  the  OH  lines  were 
obtained  also  show  far-infrared  molecular  absorption  (Kegel 
et  al.  1999;  Harvey  et  al.  1999),  the  rich  molecular  LWS  spec¬ 
trum  of  Arp  220  is  unique  for  its  high  signal-to-noise  ratio. 

The  physical  and  chemical  processes  that  account  for  the 
rich  molecular  far-infrared  spectrum  of  Arp  220  can  be  better 
understood  if  quantitative  values  of  the  column  densities  of 
the  above  species,  as  well  as  their  excitation  conditions,  are 
estimated.  The  presence  of  OH  and,  to  some  extent,  of  H20, 
may  be  indicative  of  PDR  and/or  diffuse  interstellar  cloud 
chemistry,  and  their  column  densities  potentially  give  an  es¬ 
timate  of  the  UV  field  in  the  source.  On  the  other  hand,  large 
amounts  of  H20  are  produced  in  nondissociative  shocks  (e.g., 
Cemicharo  et  al.  1999),  where  the  OH  abundance  is  also  en¬ 
hanced  (Watson  et  al.  1985).  The  OH  abundance  is  expected 
to  be  generally  higher  than  that  of  H20  in  fast  dissociative 
shocks  (Neufeld  &  Dalgamo  1989).  H20  ice  in  grain  mantles 
may  also  efficiently  return  to  the  gas  phase  through  sublima¬ 
tion  of  mantles  in  “hot  core”  regions.  But  whatever  the 
characteristics  of  the  regions  producing  the  observed  molec¬ 
ular  features,  the  lines  under  study  lie  at  wavelengths  where 
the  enormous  infrared  continuum  flux  approaches  its  maxi¬ 
mum,  and  the  molecular  excitation  of  high  dipole  species 
should  be  strongly  affected  by  absoiption  of  continuum  radi¬ 
ation.  Hence  any  reliable  estimation  of  molecular  column 
densities  require  accurate  models  for  the  dust  emission.  Un¬ 
fortunately,  ISO's  lack  of  angular  resolution  forces  us  to  rely 
on  plausibility  arguments  in  our  assumptions  about  the  regions 
where  the  different  lines  are  formed;  some  of  these  are  based 
on  general  requirements  of  excitation,  and  others  on  con¬ 
clusions  from  observations  of  galactic  sources.  The  main  goal 
of  this  work  is  thus  to  shed  light  on  the  physical  and  chemical 
processes  in  Arp  220,  based  on  detailed  model  fits  of  its 
continuum  and  far-infrared  molecular/atomic  line  absorption 
and  emission  spectrum.  We  adopt  a  distance  to  Arp  220  of 
72  Mpc  (projected  linear  scale  of  350  pc  arcsec-1;  Graham 


et  al.  1990).  In  §  2  we  present  the  ISO  LWS  observations;  in 
§  3  we  discuss  the  line  identifications;  §  4  is  a  discussion  of  the 
models  for  the  continuum  emission;  §  5  presents  the  models 
for  the  molecular  and  atomic  species;  in  §  6  we  discuss  the 
implications  of  the  radiative  transfer  models;  §  7  summarizes 
our  main  results. 

2.  OBSERVATIONS 

The  full  43-197  pm  spectrum  of  Arp  220,  obtained  with  the 
LWS  spectrometer  (Clegg  et  al.  1996)  on  board  ISO  (Kessler 
et  al.  1996)  (TDT4  number  27800202),  was  presented  by 
Fischer  et  al.  (1997,  1999).  The  grating  spectral  resolution  is 
«0.3  pm  in  the  43-93  pm  interval  (detectors  SW1-SW5)  and 
«0.6  pm  in  the  80-197  pm  interval  (detectors  LW1-LW5), 
corresponding  to  At;  >  103  km  s  1 .  The  lines  are  thus  not 
resolved  in  velocity.  The  beam  size  of  «80"  ensures  that  all 
the  continuum  and  line  emission/absorption  from  Arp  220 
(CO  size  <4";  Scoville  et  al.  1997,  hereafter  SYB97)  lie  within 
the  ISO  aperture. 

The  data  we  present  here  (co-added  TDT  numbers 
27800202,  64000916,  and  64000801)  were  reduced  using 
version  9.1  of  the  Off  Line  Processing  (OLP)  Pipeline  system, 
which  we  found  to  produce  a  higher  signal-to-noise  spectrum 
than  OLP  10.1.  However,  we  adopted  the  continuum  correc¬ 
tion  given  by  OLP  version  10.1,  which  typically  gives  ab¬ 
solute  responsivity  corrections  with  uncertainty  factors  of 
~3  times  lower  than  are  produced  by  version  9  (T.  Grundy 
2003,  private  communication).  We  performed  subsequent 
data  processing,  including  further  deglitching,  co-addition, 
and  scaling,  using  the  ISO  Spectral  Analysis  Package  (ISAP; 
Sturm  et  al.  1998).  In  order  to  obtain  a  smooth  spectrum 
throughout  the  whole  LWS  range,  the  flux  densities  given  by 
each  detector  were  corrected  by  multiplicative  scale  factors. 
Corrections  were  lower  than  10%  except  for  detector  SW1 
(43-64  pm),  for  which  the  correction  was  of  15%.  Thus  we 
attribute  a  conservative  uncertainty  of  15%  to  the  overall 
continuum  level. 

The  LWS  spectrum  of  Arp  220  is  presented  in  Figure  1 
together  with  identifications  of  the  most  prominent  lines. 
Owing  to  transient  effects,  the  fluxes  of  weak  lines  as  observed 
in  the  forward  and  reverse  scans  were  found  to  differ  signifi¬ 
cantly  in  some  wavelength  ranges.  In  these  cases,  if  a  line 
appeared  close  to  the  upper  or  lower  end  of  a  detector,  the 
reverse  or  forward  scan  was  selected  for  that  line,  respectively, 
as  has  been  found  to  minimize  the  transient  effects  (T.  Grundy 
2003,  private  communication).  Nevertheless,  the  average  of 
both  scans  was  used  throughout  most  of  the  spectrum.  In 
wavelength  regions  where  two  detectors’  responses  overlap, 
line  fluxes  were  generally  found  to  be  consistent.  The  only 
exception  was  the  H20  322-2n  line  at  90  pm,  which  showed 
in  LW1  a  flux  60%  weaker  than  in  SW5.  We  adopted  here  the 
SW5  spectrum,  but  the  flux  of  the  above  H20  line  should  be 
considered  highly  uncertain.  The  subtraction  of  a  baseline  (see 
Fig.  1)  added  additional  uncertainty  to  the  line  fluxes,  par¬ 
ticularly  in  cases  of  broad  features  presumably  composed  of 
several  lines.  With  the  exception  of  the  H20  322-2u  line,  we 
estimate  a  line  flux  uncertainty  generally  lower  than  35%. 

We  note  for  reference  that  for  the  LWS  spectral  resolu¬ 
tion  An  =  1500  km  s_1,  a  line  that  is  optically  thick  at  line 
center  with  a  line  width  of  300  km  s_l  and  negligible  line  or 


4  Target  Dedicated  Time  number,  a  unique  identifier  for  each  observation. 
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Fig.  1. — ISO  LWS  spectrum  of  Arp  220,  where  the  most  prominent  line  features  are  identified  (see  text).  The  gray  line  shows  the  adopted  baseline  (continuum 
level).  Wavelengths  in  this  and  following  figures  are  rest  wavelengths. 


continuum  re-emission,  would  show  absorption  at  line  center 
of  ~20%  in  Figure  1. 

3.  GENERAL  RESULTS 

OH  and  H20. — The  far-infrared  (FIR)  spectrum  of  Arp  220 
is  dominated  by  unresolved  OH  doublets  (hereafter  simply 
referred  to  as  lines)  and  H20  lines  in  absorption,  with  the 
exception  of  the  OH  ni/2  3/2-  1/2  emission  line  at  163.3  /mi. 
Figure  2  shows  the  level  diagram  of  OH,  ortho-H20,  and  para- 
H20  and  indicates  the  lines  detected  in  Arp  220.  Lines  with 
very  different  excitation  requirements  are  observed  throughout 
the  spectrum.  The  OH  lines  n3//2  J  =  9/2 -7/2  at  65  pm  and 
n1/2  J  =  7/2— 5/2  at  71  pm  have  lower  levels  at  290  and 
415  K  above  the  ground  state,  respectively,  whereas  the  H20 
lines  432— 32i  (59  pm)  and  422— 3 13  (58  pm)  have  lower  levels 
at  305  and  205  K.  Strong  absorption  is  also  observed  in  the 
OH  ground  state  lines  at  53,  79,  and  119  pm,  as  well  as  in  the 
H20  lowest  lying  line  at  179  pm.  This  wide  range  of  excitation 
suggests  that  several  regions  with  different  physical  conditions 
are  contributing  to  the  observed  features  (Fischer  et  al.  1999), 
and  one  of  the  goals  of  this  work  is  to  provide  a  reasonable 


estimate  of  the  nature  of  these  regions  and  their  relative  con¬ 
tributions  to  the  spectrum.  On  the  other  hand,  several  of  the 
lines  have  complex  shapes,  with  evidence  of  shoulders  sug¬ 
gestive  of  weaker  secondary  lines.  In  particular,  the  ground 
state  OH  119  pm  line  shows  a  redshifted  “shoulder,”  which 
is  detected  in  both  the  forward  and  reverse  scans,  although 
with  somewhat  different  strengths.  It  could  be  attributed  to  the 
18OH  n3/2  5/2— 3/2  line  at  120.1  pm,  although  contamination 
by  other  species  such  as  CH+  cannot  be  ruled  out.  Also,  the 
redshifted  “wing”  of  the  H20  212 — l0i  line  at  179  pm,  attrib¬ 
uted  in  Figure  1  to  the  H20  221-212  line,  could  also  be  con¬ 
taminated  by  H/sO,  CH  (see  Fig.  7),  and  H30+. 

CH  and  NH. — The  spectrum  contains  lines  from  other 
molecular  species:  like  CH  at  149  pm,  NH3  at  125,  166,  and 
170  pm  and,  very  interestingly,  strong  absorptions  at  102  and 
153.2  pm  that  have  been  identified  as  NH  in  Figure  1.  Evi¬ 
dence  for  the  latter  identifications  is  strengthened  because 
of  the  presence  of  weak  linelike  features  at  155.74  pm,  and 
marginally  at  151.53  pm,  which  would  correspond  to  the 
NH  2j—  1 !  and  2j-l0  lines,  respectively  (see  also  Fig.  3). 
Conceivably,  the  line  absorptions  at  102  and  153.2  pm  could 
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Fig.  2. — Energy  level  diagrams  of  OH  and  H2O  (ortho  and  para).  Rotational  levels  with  energies  above  the  ground  state  up  to  620  K  are  shown;  the  lines  detected 
in  Arp  220  are  indicated  with  arrows,  and  their  wavelengths  are  in  /im.  OH  A-doubling  is  ignored  because  the  A-doublets  are  not  resolved  with  the  ISO  grating 
resolution. 


be  severely  contaminated  by  other  species,  such  as  C3,  H^O. 
NH3,  and  even  OH+.  C3  has  a  strong  transition  at  153.3  //in. 
but  its  contribution  is  expected  to  be  minimal  because  of  the 
lack  of  detection  of  other  adjacent  C3  lines  (in  particular  at 
154.86  /im).  The  absorption  at  102  /im  may  be  contaminated 
by  the  H2  180  22o-ln  line  (just  at  102.0  /im),  but  since  the 
Hj,80  line  should  not  be  as  strong  as  the  corresponding  adja¬ 
cent  line  of  the  main  isotope  Hj60,  we  regard  this  identification 
as  also  unlikely.  Some  contribution  of  NH3  lines  to  the  102  /im 
feature  might  be  expected,  but  they  are  somewhat  shifted  in 
wavelength  (they  lie  between  101.5  and  101.7  /im).  Both 
features  could  be  contaminated  to  some  extent  by  OH+,  with 
strong  lines  at  101.70  and  101.92  /im  ( Nj  =  33-22  and  34-23) 
and  at  153.0  /i m  (23-l2),  but  the  strong  absorption  at  153.2  /im 


Fig.  3. — Line  assignments  in  the  vicinity  of  the  153.2  /un  feature.  The  gray 
line  shows  the  adopted  continuum  level. 


absorption  could  never  be  explained  by  OH+  alone.  Therefore, 
despite  the  possible  contamination  from  other  molecules,  NH 
is  probably  responsible  for  most,  if  not  all,  of  the  observed 

153.2  /tm  absorption  (see  also  Fig.  7).  Thus,  although  the 
definitive  assignment  to  NH  should  await  confirmation  with 
higher  spectral  resolution  observations  of  the  lowest  lying  NH 
transitions  at  ~403  GHz,  we  conclude  that  ISO  observations 
strongly  support  its  detection  in  Arp  220  and  advance  a  model 
of  the  observed  absorption  (§  5)  that  can  be  useful  to  direct 
future  observational  and  theoretical  studies.  If  confirmed,  this 
detection  is  the  first  extragalactic  detection  of  NH,  which 
has  been  previously  detected  only  in  galactic  diffuse  clouds 
through  electronic  transitions  (Meyer  &  Roth  1991;  Crawford 
&  Williams  1997)  and,  interestingly,  toward  Sgr  B2  in  the 
Galactic  center  via  the  same  transitions  detected  in  Arp  220 
(Cemicharo  et  al.  2000;  Goicoechea  et  al.  2004,  hereafter 
GRC04). 

NH3. — The  spectrum  of  Arp  220  around  125  and  170  /im, 
shown  in  Figure  4,  strongly  supports  the  identification  of 
NH3.  The  shape  of  the  165.7  /im  feature  indicates  transient 
effects,  but  the  line  is  detected  in  both  the  forward  and  reverse 
scans.  Some  H20  lines  may  contribute  to  the  observed  125 
and  127  /im  features,  but  they  are  shifted  in  wavelength  rel¬ 
ative  to  the  strongest  absorption.  This  is  the  first  extragalactic 
detection  of  infrared  NH3  lines.  With  the  detection  of  NH  and 
NH3  we  might  expect  to  detect  NH2,  but  there  is  no  evidence 
for  its  strongest  expected  lines  at  159.5,  117.8,  and  104.9  /tm. 

CO? — An  apparent  emission  line,  detected  in  both  the 
forward  and  reverse  scans,  is  present  at  173.7  /im  in  Figure  1. 
It  coincides  rather  well  with  the  expected  position  of  the  CO 
J  =  15-14  line  at  173.63  /tm.  This  identification  cannot  be 
confirmed  by  the  detection  of  other  expected  CO  lines  because 
the  CO  J  =  14-13  line  at  186.0  /im  lies  at  the  noisy  edge  of 
the  LW5  detector  and  the  CO  J  =  16-15  line  at  162.8  /mi  is 
blended  with  the  OH  163.3  /im  line.  The  higher  J  lines  are 
expected  to  be  too  weak  to  be  detectable,  given  the  observed 
strength  of  the  173.7  /im  feature. 

[O  i]  and  [C  nj. — The  ISO  spectrum  also  shows  the  [O  i] 

63.2  n m  line  in  absorption  and  the  [C  n]  157.7  /tm  line  in 
emission.  The  [O  i]  145.5  /tm  line  is  not  detected.  These  lines, 
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Fig.  4. — NH3  lines  around  125,  127,  166,  and  170  jii n.  The  gray  line  shows  the  adopted  continuum  level.  The  125  and  127  fin 1  features  could  be  partially 
contaminated  by  the  labeled  H2O  lines. 


as  observed  in  Arp  220  and  other  ULIRGs,  have  been  dis¬ 
cussed  elsewhere  (Fischer  et  al.  1997,  1999;  Luhman  et  al. 
1998,  2003).  In  §  5  we  present  a  simple  model  of  the  Arp  220 
spectrum  that  may  shed  some  light  on  the  peculiar  behavior  of 
these  lines  in  Arp  220. 

In  summary,  the  far-infrared  spectrum  of  Arp  220  shows 
molecular  lines  of  OH,  H20,  CH,  NH,  and  NH3.  The  atomic 
lines  of  [O  i]  at  63  /tm  and  [C  ii]  at  158  fim  are  also  detected. 
Lines  of  other  species,  such  as  CO,  H30+,  CH+,  and  OH+, 
could  also  contaminate  the  observed  features,  but  our  limited 
spectral  resolution  prevents  the  possibility  of  unambiguous 
detection.  Only  the  [C  n]  158  /mi  and  the  OH  3/2-1/2 
163  /on  lines  are  clearly  observed  in  emission.  Lines  from  ions 
that  would  trace  H  it  regions  and/or  an  active  galactic  nucleus 
(AGN)  are  absent. 

4.  MODELS  FOR  THE  CONTINUUM 

Figure  1  shows  that  the  continuum  peaks  around  ^40- 
50  /im.  The  bulk  of  the  continuum  from  Arp  220  is  emitted  by 
heated  dust  grains.  At  1.3  mm  wavelengths,  Sakamoto  et  al. 
(1999)  showed  that  the  continuum  arises  almost  exclusively 
from  the  nuclei,  with  an  equivalent  size  of  ~0"4.  The  non- 
themial  contribution  at  1 .3  mm  is  expected  to  be  less  than  15% 
(see  Fig.  6  of  Anantharamaiah  et  al.  2000).  On  the  other  hand, 
Soifer  et  al.  (1999,  hereafter  S99)  have  shown  that  the  two 
nuclei  also  account  for  essentially  all  the  continuum  at  25  /im. 
Combining  both  observations,  S99  proposed  two  alternative 
scenarios  to  explain  the  continuum  emission  of  Arp  220  from 
far-infrared  to  millimeter  wavelengths.  Our  models  of  the 
continuum  emission  are  entirely  based  on  these  scenarios, 
which  we  have  examined  and  refined  quantitatively  on  the 
basis  of  our  ISO  45-200  /tm  spectrum. 

Model  Sj. — In  the  first  scenario  (hereafter  Si),  it  is  assumed 
that  the  emission  from  the  nuclei  is  not  significantly  attenuated 
at  25  /tm  by  foreground  material.  We  have  simulated  the 
emission  from  the  nuclei  as  arising  from  a  single  nucleus  with 
effective  size  of  0W41.  With  an  effective  dust  temperature  of 
85  K  and  optically  thick  emission  in  the  submillimeter,  the 
requirement  that  the  fluxes  at  25  and  1300  /tm  arise  from  the 


nuclei  is  fulfilled.  However,  the  emission  at  60-100  /tm,  as 
well  as  the  total  luminosity  from  the  galaxy,  is  then  under¬ 
estimated  and  a  more  extended  region  (hereafter  ER)  must 
be  invoked  to  account  for  the  remaining  flux.  We  identify 
this  surrounding  environment  with  the  extended  emission 
observed  in  CO,  HCN,  and  HCO+.  In  this  first  scenario,  then, 
most  of  the  Arp  220  luminosity  is  produced  in  the  ER,  which 
has  been  modeled  as  a  thin  disk  by  SYB97,  and  as  a  warped 
disk  by  Eckart  &  Downes  (2001).  Significantly,  if  this  model 
is  correct,  a  spatially  extended  starburst,  responsible  for  the 
bulk  of  the  far-infrared  luminosity  is  inferred. 

Our  best  fit  to  the  continuum  using  model  Si  assumptions  is 
presented  in  Figure  5a,  with  derived  physical  parameters  listed 
in  Table  1.  In  all  models,  uniform  densities  throughout  the 
different  components  are  assumed  for  simplicity.  In  Table  1 ,  Xt 
is  the  wavelength  for  which  the  nucleus  becomes  optically 
thin  (r  =1);  owing  to  the  high  opacities  involved,  the  inferred 
emission  is  rather  insensitive  to  the  spectral  index  /?.  Thus 
we  have  given  values  of  the  physical  parameters  for  /?  =  1.5 
and  2.  The  dust  mass  has  been  derived  by  assuming  a  mass- 
opacity  coefficient  of  12  cm2  g_1  at  200  /mi  (Hildebrand 
1983).  The  parameters  that  have  been  allowed  to  vary  in  our 
models  of  the  ER  are  the  dust  temperature  Tj,  the  diameter  d 
(within  the  range  500-800  pc),  and  the  dust  opacity;  /?  is  fixed 
to  2  to  ensure  negligible  emission  at  millimeter  wavelengths 
(see  S99). 

Model  S2. — In  the  second  scenario  (hereafter  S2),  the 
emission  from  the  nuclei  at  25  /mi  is  assumed  to  be  attenuated 
by  foreground  dust  with  rabs(24.5  /im)  =  1.2,  a  value  that  was 
chosen  to  be  compatible  with  the  silicate  absorption  observed 
in  S99.  The  nuclei  account  for  the  required  flux  at  24.5  and 
1300  fim  with  Tj  «  106  K,  significantly  warmer  than  the 
~85  K  temperature  in  Si,  and  as  before,  the  emission  in  the 
submillimeter  is  optically  thick  (see  Table  1).  In  this  scenario 
as  well,  however,  the  flux  at  60-100  /mi  is  again  under¬ 
estimated,  and  an  emitting  ER  must  also  be  involved  to  ac¬ 
count  for  it.  Nevertheless,  the  luminosity  from  the  warm 
nuclei  in  S2  is  enough  to  account  for  the  observed  total  lu¬ 
minosity  from  Arp  220,  so  that  the  ER  merely  reradiates  the 
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Fig.  5. — Fits  of  the  continuum  emission  from  Arp  220  in  the  25-1300  /im 
range  for  scenarios  (a)  S!  and  ( b )  S2.  Solid  line  shows  ISO  LWS  spectrum  of 
Arp  220,  and  filled  triangles  show  the  60  and  100  jim  IRAS  fluxes  for 
comparison.  Filled  squares  show  the  fluxes  measured  by  Soifer  et  al.  (1999; 
24.5  /im),  Eales  et  al.  (1989;  450  /im),  Rigopoulou  et  al.(1996;  350,  800,  and 
1100  fim),  and  Sakamoto  et  al.  (1999;  1300  /im).  Dotted  lines  indicate  the 
computed  contributions  from  the  nuclei  and  the  ER,  whereas  the  solid  gray 
line  shows  the  expected  total  flux.  The  insert  panel  shows  the  calculated  dust 
temperatures  in  function  of  the  radial  angular  distance  for  S2  (see  text). 


emission  from  the  nuclei  and  no  extended  starburst  is  then 
needed  to  provide  the  bulk  of  the  luminosity. 

Figure  5b  shows  our  best  fit  for  S2.  The  unique  parameters 
that  have  been  allowed  to  vary  in  our  models  of  the  ER  are  the 
diameter  d  and  the  dust  opacity;  (3  is  again  fixed  to  2,  and  the 
dust  temperature  Tlt  (shown  in  the  insert  panel  of  Fig.  5b)  has 
been  computed  from  the  requirement  that  the  heating  balance 
the  cooling  throughout  the  source.  The  calculation  of  Tj  is 
carried  out  by  assuming  spherical  symmetry,  with  the  nucleus, 
the  primary  heating  source,  located  at  the  center  of  the  ER. 
The  ER  is  divided  into  a  set  of  spherical  shells  to  account  for 
the  variation  of  Tj  with  the  distance  to  the  nucleus.  Once  the 
Tj  profile  is  calculated,  the  flux  contributions  of  the  attenuated 
nucleus  and  the  ER  are  computed  separately,  and  added  up  to 
give  the  total  flux.  Despite  the  good  fit  to  the  continuum  in 
Figure  5b,  S2  implicitly  supposes  a  lack  of  spherical  symmetry 


(e.g.,  the  nuclear  disk  by  SYB97)  or  some  dumpiness,  be¬ 
cause  the  derived  radial  opacity  of  the  ER  at  24.5  /im  is  11.3, 
whereas  the  adopted  opacity  of  the  absorbing  shell  in  front  of 
the  nuclei  is  rabs(24.5  /tm)  =  1.2.  Furthermore,  we  do  not  rule 
out  the  possibility  that  the  ER  is  only  partially  responsible  for 
the  foreground  dust  absorption  of  the  nuclear  emission  by 
foreground  dust.  If  the  ER  were  concentrated  in  a  thin  disk  as 
proposed  by  SYB97,  little  dust  in  the  ER  would  be  expected 
to  lie  in  front  of  the  nuclei  and  significant  dust  absorption 
would  be  attributed  to  another  component,  “the  halo”  (see 
§  5.1).  Therefore,  our  results  for  S2,  which  assume  that  the 
nucleus  and  the  ER  are  spherically  symmetric,  should  be  con¬ 
sidered  only  approximate,  but  suggestive.  The  intrinsic  geom¬ 
etry  that  underlies  S2  departs  from  spherical  symmetry  and 
implies  that  the  total  luminosity,  which  coincides  with  the 
luminosity  of  the  nuclei,  is  lower  than  the  value  inferred  in  Si, 
where  spherical  symmetry  and  uniformity  is  assumed  for  each 
component  (Table  1). 

Table  1  shows  that  (3  =  2  yields  a  dust  mass  of  ~108  M0  for 
the  nucleus  and  that  Si  gives  a  mass  1.4  times  higher  than 
S2.  Since  SYB97  infer  a  dynamical  mass  of  (6-8)  x  109  Me 
enclosed  in  the  inner  250  pc  radius,  and  since  this  region 
contains  the  nuclei  and  most  of  the  inner  disk  (the  ER),  S2 
with  (3  =  1.5  is  favored  in  our  models,  provided  that  the  gas- 
to-dust  mass  ratio  is  not  lower  than  the  standard  value  of 
«100.  The  consistency  between  the  dynamical  mass  and  the 
mass  derived  from  the  dust  emission  indicates  that  the  inter¬ 
stellar  medium  dominates  the  dynamics  in  these  inner  regions 
of  Arp  220. 

It  is  worth  noting  that  these  models  may  be  applied  to  the 
source  as  a  whole,  as  implicitly  assumed  above,  or  alterna¬ 
tively  to  each  one  of  an  ensemble  of  Nc  smaller  clouds  of 
radius  Rc  that  do  not  spatially  overlap  along  the  line  of  sight. 
The  value  of  NCR2C  determines  the  absolute  scale,  and  the  ra¬ 
dial  opacity  and  temperature  distribution  of  each  cloud  as  a 
function  of  the  normalized  radial  coordinate,  Rc,  determine  the 
continuum  shape.  Identical  results  are  found  as  long  as  the 
above  parameters  remain  constant.  Furthermore,  both  alter¬ 
natives  give  identical  total  masses  but  differ  in  the  inferred 
mean  density,  which  scales  as  \fN~c.  Here  Nc  =  1  gives  the 
lowest  mean  density  (n(H2)),  which  is  listed  in  Table  1. 
Typical  values  of  a  few  times  104  cm~3  are  derived  for  the 
nuclei,  accounting  for  the  emission  from  molecules  with  high 
dipole  moment  such  as  HCN,  HCO+,  HC3N,  and  CN  (see 
also  SYB97).  For  the  ER  we  obtain  «(H2)  <  103  cm~3;  since 
HCN  and  HCO+  appear  to  show  extended  low-level  emission 
(Radford  et  al.  1991b),  it  is  suggested  that  the  actual  density  is 
higher  than  this  lower  limit  or  that  the  gas  is  clumpy.  The  low 
density  derived  for  the  ER  may  also  be  a  consequence  of  the 
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Physical  Parameters  Derived  from  the  Continuum  Emission 
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a  Assuming  a  single  source  instead  of  an  ensemble  of  clouds. 

b  The  temperature  is  calculated  throughout  the  ER  from  the  balance  between  heating  and  cooling  with  the  nucleus  as  the  primary  heating  source. 
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spherical  shape  attributed  to  the  ER:  if  the  mass  we  derive  for 
the  ER  were  concentrated  in  a  thin  disk  with  full  thickness  of 
32  pc  (SYB97),  the  mean  FE  density  would  be  7  x  103  cm~3. 
We  adopt  the  density  given  in  Table  1,  n( H2)  =  5.3  x  1 02  cm~3, 
in  the  models  for  molecules  and  atoms,  but  we  also  explore 
the  results  obtained  with  a  density  1  order  of  magnitude  higher 
than  the  quoted  value. 

Although  both  scenarios  Si  and  S2  reproduce  the  continuum 
emission  from  Aip  220  over  the  25-1300  /mi  interval  and 
support  the  constraints  on  the  nuclear  sizes  derived  from  the 
available  high  angular  resolution  continuum  measurements,  the 
dynamical  masses  inferred  from  CO  millimeter  line  observa¬ 
tions  favor  S2  with  (3  =  1.5  for  the  nuclei.  Moreover,  as  we 
discuss  in  §§  5.3  and  6.1,  the  observed  line  absorption/emission 
also  favors  model  S2.  We  thus  adopt  scenario  S2  for  the  detailed 
modeling  and  analysis  of  the  line  emission  and  absorption. 

4.1.  Extinction 

In  both  scenarios,  the  high  brightness  and  compactness  of 
the  nuclei  imply  extreme  continuum  optical  depth,  corre¬ 
sponding  to  Ay  ~  104  mag.  This  conclusion  is  in  strong  con¬ 
trast  with  the  much  more  moderate  values  derived  from 
infrared  and  radio  hydrogen  recombination  lines  (Genzel  et  al. 
1998;  Anantharamaiah  et  al.  2000).  The  high  extinction  de¬ 
rived  here  is  the  direct  result  of  the  measured  1.3  mm  con¬ 
tinuum  flux  from  the  nuclei,  210  mJy  (Carico  et  al.  1992; 
Sakamoto  et  al.  1999),  and  the  observed  upper  limit  of  the 
corresponding  source  size,  ~0!'4  (Sakamoto  et  al.  1999). 
These  values  imply  rj3  mm  7’(;(K)  ~  40,  which  shows  that  even 
assuming  unexpectedly  high  average  dust  temperatures  (e.g., 
Tii  =  200  K)  and  (3  =  1 ,  the  dust  emission  is  still  optically 
thick  even  at  200  /tm.  On  the  other  hand,  the  radio  recombi¬ 
nation  lines  observed  by  Anantharamaiah  et  al.  (2000)  are  not 
affected  by  dust  obscuration,  although  their  predicted  fluxes 
and  the  derived  extinction  may  be  somewhat  model  depen¬ 
dent.  These  very  different  extinction  values  may  be  under¬ 
stood  if  we  assume  that  the  observed  H  recombination  lines, 
tracing  primarily  star  formation,  are  formed  in  the  outermost 
regions  of  the  nuclei,  while  a  buried  central  energy  source, 
responsible  for  the  heating  of  dust  in  the  innermost  regions  of 
the  nuclei  and  a  significant  fraction  of  the  galactic  luminosity, 
is  weak  in  recombination  lines.  If  weak  in  recombination 
lines,  the  buried  energy  source  is  presumably  weak  in  PAH 
features  and  PDR  lines  as  well,  consistent  with  the  strong 
[C  11]  deficit  in  Arp  220  (see  also  §  5).  Dust-bounded  ionized 
regions,  in  which  most  of  the  Lyman  continuum  from  nuclear 
starbursts  or  AGN  is  absorbed  by  dust  rather  than  by  gas,  may 
explain  these  properties  of  Arp  220,  as  was  proposed  by 
Luhman  et  al.  (2003).  We  further  argue  that  the  Lyman  con¬ 
tinuum  luminosities  derived  from  recombination  lines  do  not 
empirically  rule  out  the  possibility  that  an  AGN  accounts  for 
more  than  ~50%  of  the  bolometric  luminosity  of  Arp  220, 
because  of  the  high  dispersion  of  Lbol/LLyc  values  shown  by 
both  starburst  galaxies  and  AGN,  the  range  of  ZLyc  values 
derived  from  different  tracers  in  Arp  220,  and  the  uncertain¬ 
ties  in  the  assumed  extinction  law  and  the  derived  extinction 
(Genzel  et  al.  1998).  Our  derived  /V(H2)  ~  1025  cm~2  is  high 
enough  to  obscure  a  source  of  high  5-10  keV  luminosity  from 
one  or  both  nuclei  in  Arp  220,  so  that  a  hidden  AGN  is 
allowed  despite  the  relatively  weak  X-ray  luminosity  observed 
in  Arp  220  (Clements  et  al.  2002).  Haas  et  al.  (2001)  have  also 
argued  that  a  hidden  AGN  powers  much  of  the  luminosity  of 
Arp  220  on  the  basis  of  the  observed  submillimeter  continuum 
excess  relative  to  the  7.7  /tm  PAH  flux. 


5.  MODELS  FOR  MOLECULES  AND  ATOMS 
5.1.  Comparison  with  Sgr  B 2 

The  comparison  of  the  spectrum  of  Arp  220  with  that  of 
some  well-studied  galactic  sources  provides  important  clues 
about  the  regions  where  the  observed  lines  are  formed  while 
emphasizing  the  unique  features  that  characterize  the  extra- 
galactic  source.  In  this  sense,  Sgr  B2  (component  M)  is  an 
ideal  comparison  source,  because  it  shares  common  obser¬ 
vational  properties  with  Arp  220  despite  the  obvious  differ¬ 
ences  in  spatial  scale  (and  indeed  possibly  in  nature).  Figure  6 
shows  the  continuum-normalized  spectra  of  Sgr  B2  (M) 
(kindly  provided  by  J.  R.  Goicoechea)  and  Arp  220.  Sgr  B2 
harbors  newly  bom  OB  stars,  ultracompact  H  11  regions,  and 
hot  cores,  enshrined  in  a  dusty  envelope  that  is  heated  by 
shocks  and  by  radiation  and  that  radiates  a  high  1R  luminos¬ 
ity  (see  GRC04  and  references  therein).  The  envelope  has 
the  highest  extinction  in  the  direction  of  the  N  and  M  con¬ 
densations,  with  optically  thick  emission  in  the  far-infrared  up 
to  ~200  /mi,  and  with  foreground  absorption  lines  of  OH, 
H20,  CH,  and  [O  1]  observed  with  the  ISO  LWS  grating 
(GRC04,  Fig.  6).  Fabry-Perot  observations  of  Sgr  B2  have 
also  allowed  the  detection  of  other  molecular  species  such  as 
NH,  NH3,  NH2,  HD,  H30+,  and  C3,  as  well  as  high-excitation 
lines  of  OH  and  H20  (GRC04  and  references  therein).  As 
shown  above  (§  4),  the  extinction  toward  the  nuclei  of  Arp  220 
is  also  very  high,  although  the  dust  is  significantly  warmer 
than  in  the  Sgr  B2  envelope.  The  comparison  between  both 
sources  is  meaningful  (at  least  as  a  first  approximation,  and 
from  the  point  of  view  of  the  radiative  transfer)  as  long  as 
the  nuclear  region  of  Arp  220  can  be  considered  an  ensemble 
of  continuum-thick  molecular-rich  clouds  such  as  Sgr  B2.  For 
an  ensemble  of  Sgr  B2-like  clouds,  since  both  the  continuum 
and  the  line  absorption  scale  with  the  number  of  clouds,  the 
continuum-normalized  spectrum  is  the  same  as  that  of  one 
individual  cloud.  This  result  applies  even  if  the  lines  in  the 
ensemble  are  broadened  relative  to  the  one-cloud  emission 
due  to  cloud-to-cloud  velocity  dispersion  and  rotation,  pro¬ 
vided  that  the  lines  remain  unresolved  with  the  grating  res¬ 
olution.  Thus  the  differences  between  the  two  spectra  in 
Figure  6  reveal  real  differences  in  excitation  and/or  column 
densities. 

The  high-excitation  lines  of  OH  and  H20  are  much  stronger 
in  Arp  220  than  in  Sgr  B2  (M)  (see  Fig.  6).  In  particular,  the 
OH  n3/2  J  =  9/2— 7/2  65  pm  line,  with  strong  absorption  in 
Arp  220,  is  not  detected  in  the  grating  spectrum  of  Sgr  B2,  and 
the  OH  n3/2  J  =  7/2 -5/2  84  pm  line  is  also  much  weaker  in 
Sgr  B2  (M).  This  strongly  indicates  the  presence  in  Arp  220 
of  a  high-excitation  region  with  relatively  high  OH  column 
densities.  The  high  spectral  resolution  (At;  ~  35  km  s  1 )  of 
LWS  Fabry-Perot  mode  observations  allowed  Goicoechea  & 
Cemicharo  (2002)  to  detect  high-excitation  OH  lines  in 
Sgr  B2  and  showed  that  they  are  pumped  through  absorption 
of  far-infrared  photons.  In  Arp  220  also,  these  lines  appear  to 
be  pumped  by  the  strong  infrared  radiation  flux  in  the  neigh¬ 
borhood  of  the  nuclei  (§  5.2).  Toward  Sgr  B2  (M),  Goicoechea 
&  Cemicharo  (2002)  derived  A(OH)  «  2  x  1016  cm-2,  and  we 
may  expect  significantly  higher  column  densities  toward  the 
nuclei  of  Arp  220.  The  peculiarity  of  Arp  220  is  also  revealed 
by  the  relatively  strong  absorptions  in  the  NH3  and  NH  lines. 
These  species  have  also  been  detected  toward  Sgr  B2  (M) 
with  Fabry-Perot  spectral  resolution  (Ceccarelli  et  al.  2002; 
GRC04),  but  Figure  6  indicates  much  higher  column  densities 
in  Arp  220,  at  least  toward  the  nuclei. 
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Fig.  6. — Continuum-normalized  spectra  of  Sgr  B2  (M)  and  Arp  220.  The  main  carriers  of  some  line  features  are  indicated. 


GRC04  found  that  the  OH  Ul/2  J  =  3/2- 1/2  163  /im  line, 
pumped  through  absorption  of  photons  in  the  OH  n1/2-n3/2 
J  =  3/2 -3/2  53.3  /im  line,  shows  emission  over  a  large  re¬ 
gion  associated  with  Sgr  B2.  In  Aip  220  the  line  is  strong, 
suggesting  significant  widespread  emission,  i.e.,  from  the 
ER.  It  is  also  worth  noting  that,  although  the  [C  n]  line  is 
not  detected  in  the  grating  spectrum  of  Sgr  B2  (M),  Fabry- 
Perot  observations  allowed  its  detection  (Vastel  et  al.  2002; 
GRC04),  with  a  flux  of  «1.5  x  10~17  W  cm~2  for  the  com¬ 
ponent  observed  in  emission.  This  value  is,  within  a  factor  of 
2,  similar  to  strengths  in  the  surrounding  region,  where  the 
continuum  is,  however,  much  weaker  (GRC04).  Therefore,  in 
addition  to  effects  of  self-absorption  and  absorption  of  the 
continuum  by  C+  in  foreground  excitation  clouds  (Vastel  et  al. 
2002;  GRC04),  the  low  [C  ii]/FIR  ratio  at  Sgr  B2  (M)  is  due 
to  a  strong  increase  in  the  FIR  emission  that  is  not  accom¬ 
panied  by  a  corresponding  rise  of  PDR  line  emission.  Dust- 
bounded  ionized  regions  (Luhman  et  al.  2003),  together  with 
extinction  effects  in  the  far-infrared,  might  account  for  the  lack 
of  PDR  line  emission  associated  with  this  additional  infrared 
component.  This  is  similar  to  our  hypothesis  for  Arp  220’s 
nuclei  (§  5.6).  The  [C  n]  emission  in  Arp  220  is  expected 
to  arise  from  PDRs  in  the  ER,  where  the  bulk  of  the  PAH 
emission  is  found  (Soifer  et  al.  2002).  As  shown  in  Figure  5, 
the  ER  dominates  the  observed  far-infrared  emission,  but 
since  the  main  heating  sources  are,  according  to  model  S2,  the 
nuclei,  the  ER  is  mainly  heated  via  absorption  of  infrared 
rather  than  UV  photons,  its  intrinsic  luminosity  is  relatively 
low,  and  hence  the  [C  n]  line  remains  weak  (§  5.6). 

Finally,  we  stress  that  Arp  220  and  Sgr  B2  (M)  show 
similar  strong  absorptions  in  the  ground-state  lines  of  OH  and 
H20,  as  well  as  in  the  CH  line  at  149  /im  and  in  the  [O  i]  at 
63  /im.  Fabry-Perot  observations  of  these  lines  toward  Sgr  B2 


(M)  (Cemicharo  et  al.  1997;  Stacey  et  al.  1 987;  Baluteau  et  al. 
1997;  Vastel  et  al.  2002;  GRC04),  indicate  that  most  of  the 
absorption  is  produced  at  velocities  from  —100  to  30  km  s_1 
associated  with  diffuse  low-excitation  clouds  located  along  the 
line  of  sight  to  Sgr  B2  but  not  physically  associated  with  it. 
The  similar  absorptions  found  in  Arp  220  strongly  suggest  the 
presence  of  such  a  diffuse,  absorbing  medium.  In  fact,  we  have 
found  that  the  combination  model  of  the  nuclei  and  the  ER, 
which  reasonably  reproduces  the  high-excitation  OH  and  H20 
lines,  fails  to  explain  the  strong  absorptions  in  the  lowest  lying 
OH  and  H20  lines.  We  suggest  that  in  Arp  220,  this  diffuse, 
absorbing  component  is  also  responsible  for  the  low  ratio  of 
the  CO  (2-1)  line  to  1.3  mm  continuum  emission  toward  the 
nuclei  of  Arp  220  (Sakamoto  et  al.  1999),  and  for  the  low 
brightness  of  both  the  CO  (1-0)  and  CO  (2-1)  lines  (SYB97), 
which  indicates  line-of-sight  blocking  of  the  nuclear  CO 
emission  by  low-excitation  gas  (Sakamoto  et  al.  1999).  The 
above  points  strongly  suggest  the  presence  of  an  absorbing 
component  in  front  of  the  nuclear  region  in  which  both  the 
particle  density  and  the  infrared  radiation  density  are  relatively 
low.  This  component  could  also  account  for  significant  ab¬ 
sorption  of  the  nuclear  continuum  emission  in  scenario  S2  if 
the  ER  were  a  thin  disk  (§  4).  We  refer  to  this  diffuse  com¬ 
ponent  of  Arp  220  as  “the  halo.” 

5.2.  Outline  of  the  Models 

The  dust  models  described  in  §  4  set  up  the  basis  for  the 
molecular  calculations.  These  are  carried  out  with  the  method 
described  in  Gonzalez-Alfonso  &  Cemicharo  (1997,  1999), 
which  computes  the  statistical  equilibrium  populations  of  a 
given  molecule  by  assuming  spherical  symmetry  and  line 
broadening  caused  by  microturbulence  and/or  radial  velocity 
gradients.  In  the  present  calculations  we  have  assumed,  for 


No.  1,  2004 


FAR-INFRARED  SPECTRUM  OF  ARP  220 


255 


simplicity,  pure  microturbulent  line  broadening,  but  some  tests 
showed  that  the  inclusion  of  a  radial  velocity  gradient  hardly 
modified  the  results.  Rotational  motion,  which  is  present 
around  the  nuclear  region  of  Arp  220  (e.g.,  SYB97;  Sakamoto 
et  al.  1999;  Downes  &  Solomon  1998),  is  not  included. 
Nevertheless,  the  assumption  of  spherical  symmetry  and  our 
neglect  of  steep  velocity  gradients  that  may  result  from  cloud- 
to-cloud  velocity  dispersion  may  be  considered  more  critical. 

Our  nonlocal  code  accounts  for  radiative  trapping  in  the 
molecular  lines,  collisional  excitation,  and  excitation  through 
absorption  of  photons  emitted  by  dust.  The  dust  parameters 
derived  above  for  Si  and  S2  are  used  in  the  calculations  for 
molecules.  As  in  the  case  of  the  continuum  models,  we  have 
modeled  the  nucleus  and  the  ER  separately.  This  is  required, 
in  S2,  by  the  relatively  low  dust  opacities  of  the  ER  toward 
the  nucleus  as  compared  with  the  radial  opacity  of  the  ER  (see 
§  4),  and  it  inevitably  involves  an  additional  uncertainty. 

Since  the  dust  in  the  nucleus  is  very  optically  thick 
throughout  the  ISO  wavelength  range,  model  results  are  sen¬ 
sitive  to  the  molecular  column  densities  only  in  the  outer  parts 
of  the  nucleus.  It  is  therefore  assumed  that  the  molecular 
abundances  are  different  from  zero  only  in  the  outer  regions 
of  the  nucleus,  where  the  infrared  lines  are  formed  (see 
Goicoechea  &  Cemicharo  2002  for  the  case  of  Sgr  B2).  Since 
dust  and  molecules  are  assumed  to  be  coexistent,  extinction 
effects  within  the  nucleus  are  implicitly  taken  into  account; 
they  place  important  constraints  on  the  molecular  abundances. 
We  have  adopted  a  molecular  shell  thickness  of  2  x  1018  cm, 
which  for  a  mean  n(H2)  =  4.6  x  104  cm~3  (in  S2)  corresponds 
to  Ay  ~  50  mag  and  r(50  pm)  «  0.3.  For  lines  around  50  pm, 
the  contribution  to  the  absorption  by  molecules  located  deeper 
into  the  nucleus  was  checked  in  some  tests  to  be  relatively 
weak,  owing  to  both  dust  and  molecular  line  optical  depth 
effects.  Foreground  extinction  in  S2  was  also  taken  into  ac¬ 
count.  In  the  models  for  the  ERs,  which  have  much  lower 
continuum  opacities,  we  have  assumed  that  dust  and  molecules 
have  uniform  abundance  ratio  throughout  the  whole  region. 
The  presence  of  the  central  nucleus  is  included  in  the  calcu¬ 
lation  of  the  statistical  equilibrium  populations  but  ignored  in 
the  calculation  of  the  emergent  fluxes  to  avoid  accounting  for 
it  twice. 

Toward  the  nucleus,  absorption  of  continuum  radiation 
determines  the  excitation  of  OH  and  H20;  the  radiative  rates 
are  much  higher  than  the  collisional  ones  even  for  very  high 
pressure  molecular  gas,  such  as  is  found  in  molecular  shocks 
[«(H2)  =  5  x  106  cm~3  and  7).  =  300  K;  hereafter  we  refer 
to  these  values  as  “shock  conditions”  for  brevity].  The  rate 
coefficients  of  Offer  et  al.  (1994)  and  Green  et  al.  (1993)  were 
used  to  check  the  collisional  excitation  of  OH  and  H20,  re¬ 
spectively.  If  widespread  shock  conditions  were  present,  only 
the  absorption  of  the  lowest  lying  lines  would  be  significantly 
affected.  Since  we  use  a  halo  to  match  these  lines,  we  cannot 
distinguish  between  shock  and  nonshock  conditions  (i.e.,  the 
line  ratios  are  not  sensitive  to  m(H2)  and  Tk  within  plausible 
values).  For  simplicity,  physical  parameters  for  the  halo  are 
derived  by  assuming  nonshock  conditions  for  the  nucleus. 
Concerning  the  ER,  widespread  shock  conditions  are  not  ap¬ 
plicable  because  they  would  involve  strong  emission  in  the 
H20  212 — l0i  and  3o3~2i2  lines  and  in  most  OH  lines  (hard 
to  cancel  by  any  halo).  Our  molecular  data  are  therefore 
sensitive  only  to  the  radial  molecular  column  density  N  and 
the  microturbulent  velocity  dispersion  av,  so  that  only  these 
two  computational  parameters  are  required  to  define  a  model 
of  a  given  component. 


Similar  to  the  dust  models,  models  for  molecules  may  be 
applied  to  the  source  as  a  whole,  or  alternatively,  to  each  one 
of  an  ensemble  of  smaller  clouds  that  do  not  spatially  overlap 
along  the  line  of  sight.  Besides  the  scaling  relationships  pointed 
out  above,  the  molecular  column  density  must  remain  the 
same  to  obtain  identical  results  when  varying  Nc.  On  the  other 
hand,  once  the  number  of  continuum  sources  Nc  is  fixed  (e.g., 
Nc  =  1),  the  same  line  fluxes  are  obtained  if  both  N  and  av  are 
divided  by  the  same  factor  fc,  and  the  resulting  fluxes  are  then 
multiplied  by  fc.  The  latter  reflects  the  approximate  equiva¬ 
lence  between  one  absorbing  cloud  with  column  density  N 
and  velocity  dispersion  av,  and  fc  clouds,  with  parameters  N //. 
and  (Jv/fc,  which  overlap  on  the  sky  but  not  in  the  line-of- 
sight  velocity  space.5 

Variations  in  N  and  a,  have  different  effects  on  the  line 
absorptions.  If  av  decreases  the  absorptions  are  weaker  for 
optically  thick  lines,  but  because  of  the  increase  of  line 
opacities,  the  high-energy  levels  become  more  populated  and 
the  above  weakness  is  more  pronounced  for  the  low-lying 
lines.  On  the  other  hand,  the  increase  of  N  has  little  effect  on 
very  optically  thick  lines  (most  of  them  low-lying  lines)  but 
a  larger  effect  on  those  with  moderate  opacities. 

We  have  generated  a  grid  of  models  for  the  nucleus  and  the 
ER  of  both  Si  and  S2,  by  varying  the  above  free  parameters  N 
and  av.  In  each  model,  the  molecular  shell  is  divided  into  a  set 
of  subshells  in  order  to  account  for  the  spatial  variation  of 
the  excitation  temperatures  of  the  lines.  First  we  searched  for 
the  nucleus+ER  combination  model  that  best  matches  the 
OH  and  H20  non-ground-state  lines,  with  the  same  value  of 
av  for  both  species.  As  pointed  out  above,  the  relatively  deep 
absorptions  of  the  OH  and  H20  ground-state  lines  could  not 
be  fitted  satisfactorily  by  any  model,  and  a  halo  component 
was  added  to  match  these  lines.  The  halo  was  assumed  to  be  a 
purely  absorbing  shell;  although  the  equilibrium  populations 
were  computed  in  spherical  symmetry  assuming  a  size  of  3 
times  that  of  the  ER,  limb  emission  (i.e.,  emission  for  impact 
parameters  that  do  not  cross  the  continuum  source)  was  ig¬ 
nored  in  the  calculation  of  the  fluxes  emergent  from  the  halo. 

Once  the  model  for  OH  and  H20  is  determined,  models  for 
the  other  detected  species  are  performed,  also  keeping  fixed 
the  value  of  <r„  derived  above  for  each  component.  Figure  7 
compares  the  derived  model  spectrum  and  the  observed  one, 
and  Table  2  lists  the  inferred  parameters.  The  next  sections 
are  devoted  to  explaining  the  details  of  these  calculations  and 
results. 

5.3.  OH  and  H20 

Our  best  model  fits  for  the  high-excitation  lines  of  OH  and 
H20  come  from  model  S2,  and  the  results  presented  here  focus 
on  this  scenario.  The  main  difference  between  Sj  and  S2 
consists  in  the  higher  column  densities  and/or  broader  line 
widths  required  by  Si  to  reproduce  the  lines,  owing  to  the  fact 
that  Tj  is  significantly  lower  in  S|.  Models  for  the  nucleus 
with  broad  line  widths  (a„  >  60  km  s  ';  see  below),  however, 
predict  strong  absorptions  in  some  OH  and  H20  lines,  such 
as  the  OH  n1/2  J  =  5/2— 3/2  98  fim  and  the  H20  3 13— 202 
138  /im  lines,  which  are  not  observed.  Therefore  S2,  which 
still  requires  high  column  densities,  is  a  better  fit  to  the  data. 

Our  best  model  fits  for  the  high-excitation  lines  involve 
column  densities  of  (2-6)  x  1017  cnr2  for  both  OH  and  H20, 

5  In  fact  the  models  could  be  defined  in  terms  of  the  two  independent 
variables  N /av  and  Nfc\  nevertheless  we  will  use  the  variables  N  and  av  with 

fc  =  1. 
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Fig.  7. — Continuum-subtracted  spectrum  of  Arp  220,  compared  with  the  result  of  the  model  (in  gray).  The  line  features  that  contribute  more  to  the  model  are 


identified  (see  text). 


and  av  =  50-30  km  s-1,  toward  the  nucleus.  The  model  for 
the  nucleus  reproduces  nearly  the  whole  absorption  in  the 
OH  n1/2  7/2— 5/2  and  H | /2  lf3 /2  5/2S/2  lines,  most  of  the 
OH  n3/2  9/2— 7/2,  and  significant  absorption  in  the  other  lines, 
but  the  modeled  ground-state  lines  are  far  too  weak.  It  also 
reproduces  the  frill  absorptions  in  the  H20  422— 3 13,  432-321, 
330-22i,  33 i— 220,  and  423— 3 32  lines,  and  significant  absorption 

TABLE  2 

Column  Densities  Derived  toward  Arp  220 


Species  Nucleus  ER  Halo 

OH .  2.0  x  10  2.0  x  1017  1.8  x  l(l"' 

H20 .  2.0  x  1017  3.0  x  1016  1.5  x  1016 

CH .  ...  ...  2.0  x  1015 

NH .  1.3  x  1016  ...  1.8  x  1015 

NH3 .  2.8  x  1016  ...  4.0  x  1015 

C+  .  2.0  x  1019  2.0  x  1019 

0 .  5.0  x  1019  5.0  x  1019  2.7  x  1018 


in  the  323— 212  line  and  others.  The  somewhat  low  value  de¬ 
rived  for  <j„  is  required  to  keep  the  absorptions  weak  in  some 
low-excitation  lines  that  are  marginally  or  not  detected.  Since 
low  <jv  implies  low-velocity  coverage  for  absorption  of  the 
continuum,  the  column  densities  that  are  needed  to  explain  the 
absolute  values  of  the  absorptions  in  high-excitation  lines  are 
relatively  high.  We  stress  that  this  value  of  av  must  be  inter¬ 
preted  as  a  strong  lower  limit  on  the  line  widths  that  would  be 
observed  with  high  enough  spectral  resolution,  because  rota¬ 
tion  and  cloud-to-cloud  velocity  dispersion  would  broaden 
the  observed  lines.  The  CO  kinematic  models  of  Arp  220 
by  SYB97  found  a  generic  value  of  av  =  90  km  s_1;  as  the 
authors  discuss  this  value  should  be  considered  the  joint  effect 
of  the  local  line  width  and  the  cloud-to-cloud  velocity  dis¬ 
persion  over  a  scale  of  ~  100  pc.  Our  av  is  the  local  line  width 
involved  in  the  calculation  of  opacities  and  directly  related 
to  the  molecular  excitation,  and  thus  the  kinematic  value  of 
SYB97  must  be  considered  here  an  upper  limit. 

It  is  worth  noting  that  the  above  column  densities  are  de¬ 
rived  by  forcing  the  OH  and  H20  high-excitation  lines  to  arise 
in  the  same  region.  A  slight  improvement  to  the  fit  of  the 
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H20  lines  is  obtained  with  even  lower  av,  25  km  s_1,  and 
A?(H20)  ~  101S  cm-2.  This  may  indicate  that  H20  and  OH 
toward  the  nucleus  do  not  arise  in  the  same  regions.  We  will 
adopt  in  the  following  the  nucleus  model  with  av  =  50  km 
s  ,  corresponding  to  the  spectrum  shown  in  Figure  7  and  the 
parameters  given  in  Table  2,  and  estimate  an  uncertainty  of  a 
factor  of  ~3  on  the  derived  column  densities. 

Once  the  nuclear  component  is  fitted,  a  search  for  the 
combination  nucleus+ER  that  best  accounts  for  the  remaining 
flux  in  non-ground-state  lines  is  carried  out.  For  densities  in 
the  ER,  m(H2)  <  104  cm  3  (§  4)  and  Tk  =  Td,  the  OH  and  H20 
collisional  excitation  is  found  to  be  negligible  in  comparison 
with  the  radiative  excitation.  We  have  also  explored  the  most 
plausible  situation  that  the  OH  lines  are  formed  within  the  C+ 
region  of  PDRs  (see  §  5.6):  assuming  Tk  =  300  K  (the  max¬ 
imum  allowed  by  the  collisional  rates  of  Offer  et  al.1994)  and 
«(H2)  <  104  cm  3,  radiative  rates  still  dominate  over  colli¬ 
sional  rates,  and  results  are  found  indistinguishable  from  those 
obtained  with  lower  Tk  values.  Only  densities  above  105  cm-3 
with  Tk  =  300  K  would  give  results  significantly  different  for 
the  ground-state  transitions  of  OH. 

The  ER  mainly  accounts  for  the  OH  ni/2  3/2- 1/2  163  /tm 
line,  which  is  uniquely  observed  in  emission  (but  predicted 
in  absorption  in  the  nuclei),  for  more  than  half  of  the  absorp¬ 
tion  observed  in  the  OH  n3/2  112-512  line,  for  re-emission  in 
the  n1/2  5/2 — 3/2,  and  for  significant  absorption  in  the  three 
ground-state  OH  lines.  We  estimate  for  the  ER  /V(OH)  ~ 

2  x  1017  cm  ,  with  av  =  50  km  s_1  throughout  most  of  the 
ER,  and  av  =  90  km  s  1  just  around  the  nucleus.  This  higher 
value  of  a„  was  required  to  obtain  significant  re-emission  in 
the  nondetected  OH  1 1 ]  2  5/2-312  line;  since  here  geometrical 
effects  may  be  important,  this  result  should  be  considered  with 
caution.  For  H20  we  obtain  a  significantly  lower  N(  H20)  ~ 

3  x  1016  cm~2,  giving  significant  absorption  in  the  323— 3 12, 
220— 1  n,  and  221-110  lines,  and  some  re-emission  in  the  303- 
2 12  line.  Since  the  halo  also  yields  some  absorption  in  these 
H20  lines  but  much  deeper  absorption  in  the  ground-state 
2i2— l0i  one  (see  below),  the  relative  H20  column  density  in 
these  components  is  not  well  determined. 

In  the  halo,  the  values  of  A(OH)  and  a,  were  determined 
by  fitting  the  missing  absorption  in  the  three  OH  ground- 
state  lines.  A  value  of  av  relatively  low,  15-20  km  s-1,  and 
7V(OH)  ~  2  x  1016  cm  2,  were  found  to  reasonably  fit  the 
fluxes  of  the  79  and  119  /tm  lines,  although  the  flux  of  the 
53  n m  line  is  somewhat  underestimated  (Fig.  7).  Higher 
values  of  av  would  predict  too  much  absorption  in  the  already 
saturated  119  /tm  line.  For  H20  we  find  A(H20)~1.5x 
10 16  cnr2.  Owing  to  the  strong  radiation  field  from  the  nu¬ 
cleus  and  ER,  H20  in  the  halo  is  still  significantly  excited, 
thus  yielding  also  some  absorption  in  the  220-ln,  221-110, 
and  303-212  lines. 

Despite  the  generally  satisfactory  fit  obtained  for  the  OH 
and  H20  lines,  the  shoulder  of  the  OH  119  /tm  line,  pre¬ 
sumably  produced  by  the  18OH  n3/2  5/2-312  line,  is  not  re¬ 
produced.  The  models  assume  16OH/18OH  =  500;  however, 
the  isotopic  abundance  ratio  in  Arp  220  required  to  reasonably 
fit  the  120  /tm  shoulder  is  as  low  as  16OH/18OH  ~  50.  This 
value  cannot  be  ruled  out  if  14N  is  efficiently  converted  into 
180  in  nuclear  processing  of  high-mass  stars  and  then  effi¬ 
ciently  ejected  to  the  interstellar  medium  through  stellar  winds 
and/or  supemovae  (Henkel  &  Mauersberger  1993).  Also,  there 
is  compelling  evidence  for  isotopic  ratios  of  150-200  in 
starburst  regions  of  nearby  galaxies  (Henkel  &  Mauersberger 
1993).  Because  of  the  low  spectral  resolution  of  the  spectrum 


and  the  possibility  that  the  feature  is  contaminated  by  other 
species,  we  do  not  attempt  to  place  useful  constraints  on  this 
ratio.  Nevertheless,  we  conclude  that  a  very  low  16OH/18OH 
abundance  ratio  likely  applies  to  Aip  220,  perhaps  indicating 
an  advanced  stage  starburst  (Henkel  &  Mauersberger  1993). 

The  model  of  Figure  7  predicts  an  absorbing  flux  of 
1.8  x  10~12  ergs  s_1  cnr2  for  the  n1/2-n 3/2  512-312  line  at 
34.6  /tm,  in  reasonable  agreement  with  the  observed  flux 
of  2.1  x  10  12  ergs  s  'em  2  reported  by  Skinner  et  al.  (1997). 
Concerning  the  OH-megamaser  emission  at  1 8  cm,  our  models 
are  not  suitable  to  account  for  them.  Proper  models  for  the 
maser  emission  require  inclusion  of  the  hyperfine  splitting, 
treatment  of  far-infrared  line  overlaps,  and  a  study  of  the  in¬ 
fluence  of  the  velocity  field  (e.g.,  Randell  et  al.  1995),  which 
are  not  considered  here.  The  OH-megamaser  emission  may  be 
radiatively  pumped  through  absorption  of  photons  in  the  34.6 
and  53.3  /tm  lines,  followed  by  radiative  cascade  to  lower 
levels.  Since  the  lower  level  of  the  34.6  /tm  transition  is  the 
ground-state  n3/2  J  =  3/2  level,  we  have  found  that  about 
65%  of  the  modeled  absorption  is  predicted  to  occur  in  the 
foreground  halo  rather  than  in  the  nuclei.  As  a  consequence, 
the  pump  efficiency  that  would  be  required  in  our  model  to 
explain  the  1.667  GHz  OH-megamaser  line  via  the  34.6  /tm 
line  alone  is  higher  than  the  value  of  ~1%  derived  by  Skinner 
et  al.  (1997).  Nevertheless,  given  the  uncertainty  in  the  de¬ 
rived  nuclear  OH  column  density,  the  possibility  that  some  of 
the  maser  emission  arises  in  the  most  central  regions,  and 
the  expected  additional  contribution  to  the  pumping  by  the 
53.3  /tm  line  (Burdyuzha  &  Vikulov  1990),  we  conclude  that 
our  models  are  roughly  consistent  with  the  OH-megamaser 
excitation  scheme  discussed  by  Skinner  et  al.  (1997). 

5.4.  CH 

Since  CH  is  close  to  the  Hund’s  coupling  case  b  limit  in  its 
2n  ground  state  (e.g.,  Brown  &  Evenson  1983),  we  denote  its 
levels  through  ( N ,  J),  where  N  is  the  case  b  rotational  quan¬ 
tum  number  and  J  =  N  ±  On  the  basis  of  the  results  ob¬ 
tained  toward  the  Galactic  center  (see  GRC04),  we  assume 
that  the  (2,  3/2) — (1 ,  1/2)  CH  line  we  observe  at  149.2  /tm 
arises  in  the  halo.  We  derive  /V(CH)  «2x  1015  cm~2  by  fit¬ 
ting  the  feature.  Unlike  the  case  of  Sgr  B2,  however,  the 
submillimeter  emission  from  the  nuclear  region  of  Arp  220  is 
strong  enough  to  populate  significantly  the  (1,  3/2)  level  via 
absorption  of  photons  in  the  (1,  3/2)— ( 1 ,  1/2)  line  at  560  /tm, 
so  that  our  models  for  the  halo  predict  some  contribution 
by  CH  (2,  5/2)— (1 ,  3/2)  to  the  absorption  feature  at  181  /tm 
(Fig.  7).  The  latter  is  uncertain,  however,  because  the  feature 
at  181  /tm  could  be  contaminated  by  H30+  (GRC04)  and/or 
by  stronger  absorption  of  H280,  whose  abundance  relative 
to  H'60  could  be  enhanced  relative  to  the  assumed  value 
of  1/500.  On  the  other  hand,  there  is  a  winglike  feature  at 
118.5  /im,  observed  in  both  the  forward  and  reverse  scans, 
which  could  be  caused  by  the  doublet  CH  (3,  7/2)-(2,  5/2). 
If  so,  and  since  the  excitation  of  this  line  requires  a  rela¬ 
tively  strong  radiation  field,  there  would  be  CH  in  the  nuclei 
that  would  account  for  about  one-third  of  the  absorption  at 
149.2  /tm,  and  the  CH  column  density  in  the  halo  would  be 
two-thirds  of  the  quoted  value. 

5.5.  NH  and  NH3 

In  contrast  with  OH  and  CH,  most  of  the  column  density  of 
NH  and  NH3  that  we  model  is  contained  within  the  nuclei. 
Assuming  that  the  153.2  /im  absorption  feature  is  caused 
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entirely  by  NH,  we  have  obtained  7V(NH)  ~  1016  cnr2  toward 
the  nucleus.  The  model  reproduces  the  marginal  absorption 
feature  at  76.8  /tm,  attributable  to  the  Nj  =  4j-3j'  lines. 
For  the  above  column  density,  the  22 —  1 1  and  23-l2  lines  at 
153.2  /tm  become  saturated  and  the  associated  feature  is  not 
completely  reproduced,  so  that  we  have  added  an  additional 
halo  component  with  TV(NH)  «  2  x  1015  cnr2  (Table  2). 
Nevertheless,  this  model  still  underestimates  the  absorption 
at  102  /tm,  strongly  suggesting  the  contribution  from  other 
species  (§  3). 

Within  a  given  A'-ladder,  the  excitation  of  the  NH3  non¬ 
metastable  levels  is  determined  by  absorption  of  far-infrared 
continuum  photons,  while  the  metastable  levels  in  K  =  2,  3, . . . 
are  populated  through  collisions  (see,  e.g.,  Ceccarelli  et  al. 
2002  for  the  case  of  Sgr  B2).  In  the  model  of  Figure  7  we 
have  assumed  an  average  density  of  n(H2)  =  4.6  x  104  cm  3 
and  77  =  100  K,  but  we  have  checked  that  the  model  results 
are  insensitive  to  the  adopted  77  because  of  the  blending  of 
lines  from  different  K-ladders  to  each  spectral  feature  in  our 
spectrum  (Fig.  1).  We  obtain  7V(NH3)  ~3  x  1016  cm~2  toward 
the  nucleus  to  fit  the  absorptions  at  125  and  127  /on  (Fig.  7). 
Besides  the  absorption  in  the  lines  showed  in  Figure  4,  the 
model  predicts  significant  absorption  at « 100  and  «101.6  /tm, 
caused  by  ( J,K )  =  (5 ,K)  —  (4 ,K)  lines,  which  is  still  insuffi¬ 
cient  to  account  for  the  observed  absorption  around  101.6  /tm. 
The  model  also  fails  to  explain  the  strong  absorption  at  166  /im, 
and  therefore  an  halo  component  with  /V(NH3 )  ~4  x  1 015  cm~2 
has  been  added  to  the  global  model  of  Figure  7.  Nevertheless, 
the  halo  components  of  NH  and  NH3  should  be  considered 
uncertain,  because  variations  in  the  background  continuum 
associated  with  each  component  could  in  principle  account  for 
the  missing  flux  in  the  lines. 

Finally,  we  have  explored  the  possibility  that  the  NH2  radical 
contributes  to  the  spectrum  at  some  wavelengths.  The  expected 
strongest  absorption  from  NH2  is  found  at  « 159.5  /tm,  caused 
by  the  strongest  components  of  the  3 13— 202  ortho  line  (the 
hyperfine  structure  was  neglected  in  these  calculations,  but  the 
split  of  the  levels  due  to  the  unpaired  electronic  spin  of  2  was 
taken  into  account).  At  this  wavelength,  a  marginal  absorption 
feature  may  be  attributed  to  NH2,  and  is  approximately  fitted 
with  a  model  for  the  nucleus  where  7V(NH2)  ~  1015  cm~2.  We 
have  included  this  model  of  NH2  in  Figure  7  to  show  that  the 
expected  absorption  in  other  lines,  such  asthe  322— 2n  one  at 
105  /tm,  do  not  conflict  with  the  observations,  and  we  conclude 
that  A(NH2)  <  2  x  1015  cm-2. 

5.6.  C  11  and  O  1 

A  crucial  test  of  our  model  is  whether  it  can  reproduce  the 
[C  11]  157.7  /tm  emission  and  [O  1]  63.3  /im  absorption  lines. 
Among  ULIRGs,  Arp  220  shows  one  of  the  most  extreme 
[C  11]  deficits  (Fcii/^fir  ~  2  x  10~4;  Luhman  et  al.  2003).  The 
[C  11]  line  is  formed  within  Ay  <2  mag  from  the  surfaces 
of  PDRs  (e.g.,  Wolfire  et  al.1990),  where  the  UV  field  from 
nearby  high-mass  stars,  or  from  the  average  galactic  field 
has  not  been  significantly  attenuated.  In  this  region,  photo¬ 
dissociation  maintains  most  of  the  gas  in  atomic  or  singly 
ionized  form,  but  some  radicals,  such  as  OH  and  NH,  find  their 
maximum  abundances  there  (Sternberg  &  Dalgamo  1995).  In 
particular,  OH  is  expected  to  be  an  excellent  molecular  tracer 
of  PDRs’  surfaces,  given  that  its  abundance  is  rather  low  in 
UV-shielded  quiescent  molecular  clouds.  Its  abundance  rela¬ 
tive  to  H  nuclei  within  the  C+  region  of  dense  PDRs  is  expected 
to  approach  the  value  of  ~3  x  10  6  (Sternberg  &  Dalgamo 
1995).  In  fact,  Goicoechea  &  Cernicharo  (2002)  have  found  an 


abundance  of  ss2  x  lO-6  in  Sgr  B2,  and  it  could  be  as  high  as 
5  x  lO-6  around  the  galactic  center  (Genzel  et  al.  1985).  We 
have  estimated  W(OH)  ~(l-3)  x  lO-6  toward  the  nucleus  of 
Arp  220. 

On  the  above  grounds,  and  adopting  a  gas  phase  carbon 
abundance  of  1.4  x  lO"4  (Savage  &  Sembach  1996),  we  as¬ 
sumed  N(C+)/N(OH)  =  100  in  Arp  220  (Table  2),  and  we 
computed  the  expected  [C  11]  line  emission  by  assuming  excita¬ 
tion  through  collisions  with  atomic  H  (Tielens  &  Hollenbach 
1985).  The  collisional  rates  were  taken  from  Launay  &  Rouefif 
(1977a).  The  H  densities  were  assumed  to  be  twice  the  H2 
densities  derived  from  the  continuum  models,  i.e.,  n(H)  = 
9.2  xlO4  cm-3  in  the  nucleus  and  «(H)  =  1.06  x  103  cm-3 
in  the  ER.  Since  the  critical  density  is  3  x  103  cm-3  (Kaufman 
et  al.  1999),  results  for  the  nucleus  are  not  critically  dependent 
on  the  assumed  density.  They  are  also  insensitive  to  the  as¬ 
sumed  temperature  as  long  as  it  is  higher  than  ~  1 00  K 
(  Wolfire  et  al.  1990).  Thus  our  results  are  sensitive  only  to  the 
assumed  density  in  the  ER  and  to  the  assumed  C+  abundance. 

The  result  of  this  calculation  has  also  been  added  to  the 
overall  model  of  Figure  7  (77  =  500  K  has  been  assumed). 
The  model  overestimates  the  observed  [C  11]  emission  by  only 
24%.  The  contribution  from  the  nucleus  is  only  13%  of  that 
from  the  ER  because  of  the  low  volume  of  the  nuclear  emit¬ 
ting  region.  The  expected  line  flux  from  the  nucleus  could 
be  lower  if  absorption  of  the  underlying  continuum  by  low- 
excitation  C+  in  the  halo,  ignored  in  this  model,  occurs  as 
observed  in  Sgr  B2  (M)  (GRC04).  The  bulk  of  the  line  emis¬ 
sion  arises  from  the  ER.  If  the  density  of  the  ER  were  1  order 
of  magnitude  higher  than  assumed  (§  4),  the  modeled  [C  11] 
emission  would  be  a  factor  of  ss2  stronger  than  in  Figure  7. 
The  situation  in  Arp  220  resembles  what  is  found  in  Sgr  B2 
(GRC04),  where  the  line  is  emitted  mainly  from  an  extended 
region  around  condensations  N  and  M,  while  the  strong  FIR 
source  itself  is  not  associated  with  corresponding  observable 
[C  11]  line  emission. 

The  [O  1]  63.3  /tm  line  has  been  modeled  by  assuming 
N( O  l)//Y(OH )  =  250  (Table  2);  although  the  oxygen  abun¬ 
dance  is  expected  to  be  twice  that  of  C+  in  the  atomic  region,  it 
is  expected  that  further  atomic  oxygen  exists  deeper  into  the 
clouds  (Sternberg  &  Dalgamo  1995).  For  this  reason,  a  wide 
range  of  excitation  temperatures  is  expected  for  the  [O  1]  line. 
We  have  just  fitted  a  single  “effective”  kinetic  temperature 
and  assumed  also  collisions  with  atomic  H.  The  same  densities 
as  assumed  above  for  C+  excitation  are  used  for  O  1,  and  the 
collisional  rates  are  taken  from  Launay  &  Rouefif  (1977b). 

Our  calculations  show  that  absorption  in  the  [O  1]  63.3  /tm 
line  is  obtained,  both  toward  the  nucleus  and  the  ER,  with  an 
effective  77  =  160  K,  but  the  line  is  still  too  weak  to  account 
for  the  observed  feature.  Given  the  high  O  1  abundances 
that  are  expected  in  diffuse  clouds  (Baluteau  et  al.  1997),  we 
have  added  to  the  model  a  halo  component  with  N( O  I)  = 
3  x  1018  cnr2  (Table  2).  This  model  accounts  for  the  observed 
absorption  at  63.3  /on  (Fig.  7). 

On  the  basis  of  the  low  extinction  derived  from  infrared 
and  radio  H  recombination  lines  and  the  high  optical  depth 
derived  from  our  dust  models  (§  4.1),  and  the  assumption  that 
[C  11]  emission  is  expected  to  suffer  extinction  similar  to  that 
of  the  recombination  line  emission,  we  favor  the  explanation 
proposed  by  Luhman  et  al.  (2003;  see  also  §  5.1)  that  non- 
PDR  far-infrared  emission  is  responsible  for  the  [C  11]  deficit 
in  Arp  220.  Nevertheless,  our  continuum  models  by  them¬ 
selves  (§  4)  cannot  rule  out  the  possible  role  of  far-infrared 
extinction  on  the  measured  line  fluxes:  the  derived  high  nuclear 
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far-infrared  opacities  indicate  that  only  the  outer  regions  of  the 
nuclei,  where  the  OH  and  H2O  lines  are  formed,  and  the  ER  are 
able  to  contribute  to  the  [C  11]  line  emission.  Our  model  fits 
cannot  discern  whether  the  intrinsic  nuclear  [C  11]  emission  is 
negligible  or  rather  is  obscured  by  dust.  In  either  case,  we  have 
shown  that  the  [C  11]  line  is  well  reproduced  by  assuming  that 
C+  and  OH  are  coexistent. 

6.  DISCUSSION 

6.1.  The  Extended  Region  (“ER”) 

Our  models  support  the  widespread  presence  of  PDRs  in 
the  ER.  Both  the  high  OH  and  C  11  column  densities  indicate 
that  the  UV  field  from  newly  formed  stars  have  a  profound 
effect  on  the  chemistry  in  the  ER.  Significant  contributions 
from  shocks  can  be  neglected,  as  was  pointed  out  in  §  5.2. 
The  H20-to-OH  abundance  ratio  is  significantly  lower  than 
1,  probably  indicating  enhanced  H20  photodissociation.  As 
pointed  out  in  §  4,  the  density  is  uncertain  but  is  thought  to  be 
in  the  range  n( H2)  =  5  x  102  —  7  x  103  cm~3. 

Our  models  also  indicate  that  the  bulk  of  the  [C  11]  line 
emission  arises  in  the  ER.  It  is  therefore  likely  that  star  for¬ 
mation  is  responsible  for  this  emission.  This  result  is  strongly 
supported  by  the  observations  of  Soifer  et  al.  (2002),  who 
found  that  the  PAH  emission  is  also  spatially  extended.  As¬ 
suming  a  “normal”  [C  ii]/FIR  ratio,  i.e.,  F[C  n] /^fir  =  5  x 
1 0  3  (Stacey  et  al.  1991),  the  expected  intrinsic  FIR  emission 
from  the  ER  is  ^3  x  1010  L0,  i.e.,  ~3%  of  the  total  galactic 
luminosity.  This  estimate  supports  the  scenario  S2  that  has 
been  used  to  model  the  line  emission.  According  to  our 
models,  the  total  FIR  luminosity  (due  to  absorption  and  re¬ 
emission  of  nuclear  infrared  radiation)  from  the  ER  is  much 
higher  than  the  intrinsic  (PDR)  emission.  Our  estimate  of  the 
intrinsic  ER  luminosity  is  somewhat  lower  than  that  derived 
by  Soifer  et  al.  (2002)  and  Spoon  et  al.  (2004),  who  found  a 
lower  limit  of  ~7  x  1010  ZQ,  but  confirms  qualitatively  their 
results.  The  starburst  luminosity  in  the  ER  seems  to  be  similar 
to  that  of  moderately  bright  infrared  galaxies,  such  as  NGC 
253  (Radovich  et  al.  2001). 

The  warm  mass  in  the  ER  derived  from  our  OH  models  is 
~  I 09  Mq,  assuming  X(OH )  =  1 0  3  relative  to  H  nuclei  and 
the  warm  H2  arising  from  the  same  volume  as  the  OH.  This 
estimate  of  mass  may  be  considered  an  upper  limit.  From  the 
H2  rotational  lines  detected  by  ISO,  Sturm  et  al.  (1996)  esti¬ 
mated  ~3  x  109  Mq  of  warm  gas.  Given  the  uncertain  cor¬ 
rection  by  extinction,  it  is  likely  that  an  important  fraction  of 
the  H2  emission  arises  in  the  ER. 

6.2.  The  Halo 

The  models  also  indicate  the  presence  of  a  halo.  Despite 
the  uncertainties  in  the  column  densities  of  this  component, 
they  are  typical  of  those  found  toward  our  galactic  nucleus: 
<V(H20,  OH)  ~2x  1016  cm~2  has  also  been  derived  in  the 
diffuse  medium  toward  Sgr  B2  (Cernicharo  et  al.  1997; 
Neufeld  et  al.  2003;  Goicoechea  &  Cernicharo  2002),  and  the 
derived  TV(CH)  «2x  1015  cm~2  is  also  similar  to  that  found 
toward  SgrB2  by  GRC04  and  Stacey  et  al.  (1987).  The  column 
densities  in  the  halo  derived  for  NH  and  NH3  are  uncertain 
because  they  are  based  on  single  lines,  but  they  could  also 
exist  in  a  population  of  molecular  clouds  located  far  away 
from  the  nuclear  region.  Assuming  X(H20)  =  10~6  in  the 
halo  (Neufeld  et  al.  2000),  tV(H2)  ~2x  1022  cm~2  is  obtained 
and  the  associated  continuum  opacity  is  r2oo^m ~8x  10~3.  If 
we  further  assume  a  spectral  index  (3  =  2 ,  the  derived  dust 


opacity  at  25  /mi  is  ~0.5,  comparable  to  the  value  of  1.2 
in  scenario  S2  (§  4).  Therefore,  significant  absorption  of  the 
nuclear  continuum  emission  is  attributable  to  the  halo. 

6.3.  The  Nuclei 

The  similarity  of  the  OH  column  densities  in  the  ER  and 
the  nucleus  may  suggest  that,  at  least  to  some  extent,  we  are 
observing  the  same  widespread  OH  component,  and  that  OH 
is  more  excited  toward  the  nucleus  because  of  the  underlying 
stronger  infrared  continuum  in  that  direction.  On  this  ground  a 
PDR  origin  of  the  observed  OH  would  be  favored.  However, 
the  inferred  high  H20  column  density  would  be  difficult  to 
explain  in  this  context.  Although  H20  is  expected  to  form 
efficiently  in  UV-shielded  regions  of  dense  PDRs  (Sternberg 
&  Dalgamo  1995),  with  total  column  densities  similar  to  those 
of  OH,  there  seems  to  be  no  clear  correlation  between  MOH ) 
and  ,V(H20).  In  fact,  the  H20  column  density  in  the  ER  is 
significantly  lower  than  that  of  OH.  Also,  detection  of  NH3 
and,  above  all,  of  NH,  seem  to  indicate  an  additional  nuclear 
component. 

OH  could  also  arise  in  C-shocks  (Watson  et  al.  1985),  but  it 
is  unlikely  that  they  dominate  the  OH  absorptions  because 
;Y(H20)  would  then  be  at  least  1  order  of  magnitude  higher 
than  A(OH).  However,  an  interesting  possibility  is  that  those 
C-shocks,  or  alternatively  hot  core  regions,  are  combined  with 
PDRs,  i.e.,  H20  produced  there  is  subject  to  a  strong  UV  field 
and  thus  to  photodissociation.  This  process  may  be  enhanced 
if  there  is  hot  gas  emitting  X-rays,  such  as  in  supernova 
remnants  where  the  OH  abundance  is  expected  to  be  >  1 0  9 
(Wardle  1999),  or  from  a  nuclear  AGN.  In  the  nuclei  of 
Arp  220,  several  compact  18  cm  continuum  sources  indicate 
the  presence  of  high  luminous  supemovae  (Smith  et  al.  1998); 
however,  the  extended  and  external  OH  required  to  explain 
the  infrared  data  suggests  a  more  widespread  component.  The 
diffuse  OH  megamaser  emission  found  in  Arp  220  should  be 
related  to  it  (Lonsdale  et  al.  1998).  Soft-extended  and  hard- 
compact  X-ray  emission,  detected  around  and  from  the  nuclei 
(Clements  et  al.  2002),  could  be  responsible  for  photodis¬ 
sociation  of  H20  produced  in  shocks  and  hot  cores,  thus 
enhancing  the  OH  abundance.  In  particular,  hot  cores  are  ex¬ 
pected  to  exist  widely  in  the  nuclei,  given  the  high  dust  tem¬ 
peratures  and  densities  found  there.  The  presence  of  ./-shocks, 
where  ,V(OH )  is  expected  to  be  higher  than  V(H20)  except  for 
high  enough  preshock  densities  (Neufeld  &  Dalgamo  1989), 
cannot  be  disregarded. 

The  high  column  densities  obtained  for  NH  and  NH3 
seem  to  indicate  that  standard  gas-phase  PDR  chemistry  alone 
is  not  able  to  explain  the  full  molecular  data  in  Arp  220. 
Sternberg  &  Dalgamo  (1995)  predicted  for  a  PDR  an  NH  col¬ 
umn  density  more  than  2  orders  of  magnitude  lower  than  that 
of  OH,  whereas  we  estimate  2V(OH)/iV(NH)  ~  20  in  Arp  220. 
The  enhancement  of  NH  relative  to  OH  may  be  more  than  1 
order  of  magnitude.  In  Sgr  B2,  GRC04  also  found  a  somewhat 
high  A(NH)  relative  to  OH,  i.e.,  A(OH)/A(NH)  ~  30-100.  It 
is  interesting  that  the  high  NH  abundance  in  diffuse  clouds  is 
a  factor  of  ~40  higher  than  predicted  by  gas-phase  chemical 
models  (Meyer  &  Roth  1991).  The  latter  has  been  used  by 
Wagenblast  et  al.  (1993)  and  Crawford  &  Williams  (1997) 
to  argue  for  grain-surface  production  of  NH.  In  principle, 
this  process  could  also  help  to  enhance  the  NH  abundance  in 
Arp  220,  because  the  dust  in  the  nuclei  has  been  found  to  be 
warm  so  that  grain  mantles  could  efficiently  evaporate,  and 
also  because  in  an  environment  with  enhanced  cosmic  rays, 
as  expected  from  a  starburst  and  the  consequent  high  rates  of 
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supemovae,  or  an  AGN,  the  release  of  grain  mantles  to  the  gas 
phase  by  means  of  sputtering  should  be  also  enhanced.  How¬ 
ever,  hydrogenation  of  NH  should  in  principle  continue  until 
saturation,  because  indeed  we  observe  NH;  with  a  column 
density  twice  of  NH.  The  issue  that  now  arises  is  that,  if  hy¬ 
drogenation  generally  completes,  then  a  very  low  NH-to-NH; 
abundance  ratio  would  be  expected,  and  if  it  does  not  complete, 
like  in  the  models  of  Wagenblast  et  al.  (1993),  the  scenario 
fails  because  of  the  low  relative  abundance  found  for  NHi.  In 
Sgr  B2,  for  example,  NH:NH2:NH3  =  1:10:100  (GRC04), 
and  NH2  is  found  to  be  fairly  abundant  (van  Dishoeck  et  al. 
1993). 

One  possible  solution  for  the  low  NH2  abundance  is  that  the 
already  invoked  enhancement  of  cosmic  rays  deeply  affects  the 
ion-molecule  gas-phase  chemistry.  Fedennan  et  al.  (1996)  have 
shown  that  the  high  NH  abundance  found  in  some  diffuse 
galactic  environments  could  be  explained  through  cosmic-ray 
ionization  of  atomic  nitrogen,  followed  by  hydrogen  abstrac¬ 
tion  reactions  that  form  NH+  and  dissociative  recombination 
that  yields  NH.  The  slight  endothermicity  of  N+  +  H2  — » 
NH+  +  H  (Millar  et  al.1997)  is  not  a  problem  here,  given  the 
high  dust  temperatures  in  the  nuclei.  If  in  Arp  220  the  cosmic- 
ray  ionization  of  N  is  enhanced,  the  above  scheme  could  give 
rise  to  high  NH  abundances.  Furthermore,  both  NH2  and  NH; 
would  be  much  less  abundant  if  photodissociation  is  impor¬ 
tant  in  those  regions.  NH3  would  be  formed  primarily  in  grain 
mantles  through  nearly  complete  hydrogenation,  thus  again 
keeping  the  NH2  abundance  low,  and  released  to  the  gas  phase 
in  widespread  hot  core  regions  relatively  shielded  from  UV 
fields.  H20  could  also  follow  this  last  process. 

A  chemistry  deeply  influenced  by  ion-neutral  reactions  have 
been  also  invoked  by  Aalto  et  al.  (2002)  to  explain  the  high 
emission  from  HNC  relative  to  HCN  in  Arp  220  and  other 
luminous  infrared  galaxies.  Furthermore,  Aalto  et  al.  (2002) 
also  found  strong  subthermal  CN  emission,  indicative  of  gas 
at  moderate  densities  and  irradiated  by  UV  fields  (Rodriguez- 
Franco  et  all 998).  In  the  simplest  scenario,  the  molecular 
content  of  Arp  220  may  thus  be  interpreted  in  terms  of  hot 
cores  submitted  to  strong  UV  and  X-ray  fields,  where  enhanced 
evaporation  of  grain  mantles  and  ion-molecule  chemistry  in¬ 
duced  by  cosmic-ray  ionization  are  also  deeply  affecting  the 
relative  molecular  abundances. 

7.  SUMMARY 

Because  of  its  high  signal-to-noise  ratio,  the  rich  far- 
infrared  spectrum  of  Arp  220  provides  a  template  for  under¬ 
standing  the  dusty  interstellar  medium  of  ULIRGs.  We  have 
analyzed  the  ISO  LWS  spectrum  of  Arp  220  using  radiative 
transfer  models  applied  to  both  the  continuum  and  line  emis¬ 
sion.  Our  main  results  are  the  following: 

1.  The  continuum  emission  from  25  to  1300  /im  is  well 
reproduced  with  a  two-component  model:  (a)  the  nuclei,  with 
effective  size  of  0"4  and  dust  temperature  of  106  K,  which 
accounts  for  essentially  the  whole  flux  at  25  /tm  and  at  milli¬ 


meter  and  submillimeter  wavelengths,  and  ( b )  an  extended 
region  (ER),  whose  effective  size  is  2"  and  which  dominates 
the  continuum  emission  from  60  to  250  /t m . 

2.  The  extinction  toward  the  nuclei  is  very  high  (Ay  ~ 
104  mag);  the  dust  in  the  ER  is  heated  through  absorption  of 
radiation  emanating  from  the  nuclei. 

3.  The  spectrum  of  Arp  220  shows  molecular  lines  of 
OH,  H20,  CH,  NH,  and  NH3,  as  well  as  the  atomic  [O  i]  63  //m 
line  in  absorption  and  the  [C  n]  158  /tm  line  in  emission.  The 
outermost  regions  of  the  nuclei,  along  with  the  ER,  are  traced 
by  the  lines  observed  in  the  far-infrared. 

4.  The  high-excitation  lines  of  OH  and  H20  are  pumped 
through  absorption  of  photons  emitted  by  dust.  Column  den¬ 
sities  ofV(OH)  =  (2-6)  x  1017  cm-2  and  V(H20)  =  (2-10)  x 
1017  cm-2  are  derived  toward  the  nuclei.  In  the  ER,  AT  OH )  ~ 

2  x  1017  otT2  and  V(  H20)  ~  3  x  1016  cm-2.  We  found  it  nec¬ 
essary  to  invoke  a  third  component,  or  halo,  to  match  the 
low-lying  lines  of  OH  and  H20;  this  halo  has  column  densities 
that  are  similar  to  those  found  toward  the  Galactic  center 
[N(OH,  H20)  ~  1.5  x  1016  cm-2]. 

5.  The  CH  line  detected  in  the  far- infrared  spectrum  of 
Arp  220  is  assumed  to  arise  from  the  halo,  and  the  inferred 
column  density  is  V(CH)  ~2x  1015  cm-2.  This  value  is  also 
similar  to  that  found  toward  the  Galactic  center. 

6.  Models  for  NH  and  NH3  indicate  high  column  densities 
toward  the  nuclei,  V(NH)  ~  1.5  x  1016  cirr2  and  V(NH3)  ~ 

3  x  1016  cuT2.  The  upper  limit  found  for  the  column  density  of 
NH2  is  much  lower,  V(NH2)  <  2  x  1015  cm-2. 

7.  The  [C  n]  158  //m  line  strength  is  approximately  repro¬ 
duced  by  assuming  that  C+  is  100  times  more  abundant  than 
OH.  Our  models  predict  that  the  line  arises  mainly  from  the 
ER,  and  that  non-PDR  far-infrared  emission,  with  possible 
extinction  effects,  is  mostly  responsible  for  the  observed  [C  n] 
deficit  in  Arp  220.  The  [O  i]  63  /im  line  is  also  matched  with  an 
abundance  of  250  relative  to  OH  and  absorption  toward  the 
nuclei,  the  ER,  and  the  halo. 

8.  PDR  molecular  chemistry  plays  a  key  role  in  the  ER  and 
contributes  to  the  elevated  OH  abundance  at  the  expense  of 
H20.  Toward  the  nuclei,  however,  important  contributions 
from  hot  cores,  and  possibly  from  shocks,  is  most  plausible. 
The  nitrogen  chemistry,  and  in  particular  the  high  NH  abun¬ 
dance,  seems  to  be  strongly  influenced  by  ion-neutral  reactions 
triggered  by  cosmic-ray  ionization. 
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